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Abstract. For tri-diagonal matrices arising in the simplified Jaynes- 
Cummings model, we give an asymptotics of the eigenvalues, prove 
a trace formula and show that the Spectral Riemann Surface is ir- 
reducible. 
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1. Introduction 

We consider one-sided tri-diagonal matrices of the form L + zB, 
where 
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and study their spectra in the case where the diagonal matrix majorities 
the off-diagonal one in the sense of the following condition (or some 
version of it) 
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There is a vast literature (see [T3J EH El I2H] and the bibliography 
therein) devoted to a broad range of questions on these matrices and 
the corresponding operators in £ 2 (N). We will be concerned with the 
following three questions. 



1. Spectra Sp(L + zB). Of course, Sp(L) 
Le k = q k e k , k = l,2, 



{q k , k = 1,2, . . .} and 
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where {e^}^ is the canonical orthonormal basis in £ 2 (N). Under the 
condition (jl.2j) the spectrum Sp(L + zB) is discrete as well (see, e.g., 
Lemma 8 in [S] , or jT21 E2] ) , and 

Sp(L + zB) = {E n (z)}?, 

where, for each n, E n (z) is an analytic function at least for small \z\, 
i.e., in the disk \z\ < R n for some R„ > 0. 
(l.A) How large could R n be chosen? 

Let us mention that in the case of Mathieu operator H. Volkmer [37] 
proved that R n x n 2 (see further discussion in Section 7.1-7.3). 

(l.B) What is the asymptotic behavior of E n (z) if z is bounded, say 
\z\ < R, and n — > oo? 

2. Under the conditions (jl.2j) and some further assumptions on the 
sequences q, b, c one can introduce the regularized trace 

oo 

tr(z) = ^T(E n {z)-q n ) 

n=l 

as an entire function - see Definition in Section 5.4. 
Can we evaluate it in specific examples? 

3. Spectral Riemann Surface of the pair (L,B) G (jl.lj) . (jl.2|) is de- 
fined as 

G = {(X,z)eC 2 : (L + zB)f = Xf, fe i 2 (N), f ^ 0}. 

F. W. Schafke proved that in the case of the Mathieu equation 

— y" + z(cos 2x)y = \y, i.e., L = —(d/dx) 2 , By = (cos 2x)y, 

the Spectral Riemann Surface is irreducible [20], PP- 88-89; see also 
|4*Uj . We use Schafke's scheme to prove that the Spectral Riemann 
Surface G is irreducible in the case of the simplified Jaynes-Cummings 
model (Theorem 3). 

We focus our attention on special tri-diagonal matrices which are 
motivated by the analysis of second order differential operators in the 
framework of Fourier method. 

Example 1. Let 

(1.3) q k = k 2 } b k = c k = k a , 0<a<2. 
If 

(1.4) a = 0, 
we have the Mathieu matrices, and if 

(1.5) a = 1/2, 
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we have the simplified Jaynes-Cumming matrices that have been con- 
sidered by A. Boutet-de-Monvel, S. Naboko and L. Silva jl]. 

Example 2. More general q, 

q k = F, b k = c k = k a , 7 > a + 1/2. 

The case 7 = 1, a = 1/2 comes from the Jaynes-Cumming model (see 
E. Tur |IIH USD- 
Example 3. The Whittaker-Hill matrices (see [THJ , Ch.7, and [S]) 

(1.6) q k = k 2 or (2k + l) 2 , b k = t - k, c k = t + k, t>0 fixed. 

We do not provide details about the Fourier method or the gauge trans- 
form which lead us from the differential operator 

— y" + (a cos 2x + b cos 4x)y 

to the matrices JHU with ffTHj). See jTH EH CHI U ■ In Section 7.1, 
Propositions HH1 and Unj we use results about differential operators (23 
EE1 EH| to find asymptotics of the radius of analyticity R n in the case 
of matrices (|1.6j) . 

The matrices (|1.3|) - (|1.5|) and (|2.1j) . (|2.2j) is the main object of interest 
in this paper. Now we spotlight some of its results. Below E n (z) means 
the n-th eigenvalue of L + zB. 

Theorem 1. Suppose \2. 1\) and \2.2}) with < a < 1/2 hold, and 
limfc b k c k k~ x = £ exists for a = 1/2. Then, for a G [0, 1/2], the regular- 
ized trace tr(a, z) is well-defined entire function, and 

(1.7) tr(a, A = f>„ W -„>) = { ° 2 I <f ■ 

See further comments in Section 7.6, Proposition 1231 

Theorem 2. Suppose that M.ty) holds and a G [0, 2/3]. For each R > 
i/iere is > s^c/i that for n > the eigenvalues E n (z), \z\ < R, 
are well defined and 
(1.8) 

E n (z) = n 2 + z 2 ( i^S + + (1 ~ 2a)( ^: 14a + 3) Vo (n«^M) 

See Theorem 1111 in Section 4.4 also. (For a = 1/2 similar formula 
was given in 4 J but it was not correct). 

Theorem 3. In the case / TO]) with a G [0,0.085], or a G [(2 - 
V / 2)/4, 1/2] t/ie Spectral Riemann Surface 

G = {(A, z)gC 2 : AG Sp(L + zB)}, 
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is irreducible. 

See further comments in Section 7.5, Proposition 1221 In the case of 
anharmonic oscillator 

Ly = —y" + x 4 y, By = x 2 y, x G M 

a question about structure of SRS and its branching points has been 
raised and solved (!) by C. Bender and T. Wu [T]; see also [2E1 EH EHl 

ESI EEl EZ! ■ 

The case of Mathieu-Hill operators could be deduced to Example 
(HH+CEH; it has a longer history (see [13 [23 12 13 E3 E3 E3 HD|). 
Some observations about Whittaker-Hill operators could be found in 
[5], Section 5.4. 

4. In the course of proving Theorems ^[3] we observe a series of facts 
and inequalities about the eigenvalues of the operators L + zB which 
could be of some interest by themselves. We discuss them in detail in 
related sections of the paper or in Section 7. 

2. Localization of the spectra 

1. Well-known methods of Perturbation Theory give information 
about the spectra Sp(L + zB) if L, B 6 (jl.lj) . (|1.2|) . For a while, let us 
assume that the sequences q, b, c satisfy the conditions 

(2.1) q k = k 2 - 



(2.2) |&fc|, \c k \ < Mk a , 0<a<2. 
For each n 6 N we set 

(2.3) A n = {z e C : |s| < R n }, R n = n 1_a /(8M). 

Proposition 4. Under the conditions \2. 1)) and the spectrum of 
the operator L + zB is discrete, and for each n and z G A n there is 
exactly one eigenvalue E n (z) in the strip 

H n = {X G C : n 2 - n < Re X < n 2 + n}. 

Moreover, the function E n (z) is analytic in A n , 

(2.4) E n (0) = n 2 , 
and 

(2.5) \E n (z) - n 2 \ < n if \z\ < R n , 
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Proof. The resolvent-operator 
(2.6) 

R x = (A - L - zB)- 1 = R° x (1 - zBR° x ) 1 , where R° x = (A - L)~\ 
is well defined if 

X£Sp(L) and |*| • \\BR° X \\ < 1. 
Let K n be the open disk with center n 2 and radius n, i.e., 

(2.7) K n = {XeC: \X-n 2 \<n}. 

By (J2.8|) (see Lemma El below) we have |z| • < 1 for \z\ < R n 

and A G \ thus 

Sp{L + zB) n (H n \ K n ) = 0. 

If z = 0, then Sp(L) = {k 2 : G N}, so n 2 is the only eigenvalue inside 
the circle dK n . It is simple, and for each z G A n the operator L + zB 
has exactly one simple eigenvalue E n (z) G if„ because 

dim f — / (A — L - zB)- l d\ J = 1. 

Moreover, it is well-known that simple eigenvalues depend analytically 
on the perturbation parameter (e.g., see ^3]), and therefore, for each 
n, E n (z) is an analytic function on A n . This completes the proof of 
Proposition I3J □ 

2. The next lemma gives the estimate of the norm 

Lemma 5. Under the assumptions \2.1\) and \2.2\) . if A = x + it G 

H n \ K n then 

(2.8) \\BR\\\ < 2M max(2, 2 a )n a ~ l , Vt G K, 

(2.9) ll^ll < 2M max(2,2 a )n7|t|, t/ n < \t\ < n 2 , 

(2.10) ll-B^II < 4M2°|t| (a - 2)/2 , i/ |t|>n 2 . 

Proof. Since = {1/(A — A; 2 )} is a diagonal operator, while B is an 
off-diagonal one, the norm ll-B-R^II does not exceed, in view of (|2.2j) . 

2Mfc Q 



(2.11) < sup 1 ' * 11 < 



k 





1 + 


Ck- 






A - 


k 2 \ 





sup 



A -PI' 



For every t G R, if k < n, then |A — k 2 \ > n — 1 > n/2, and therefore, 
k a /\X — k 2 \ < 2n a ~ 1 . For k = n we have n°7|A — n 2 | < n a ~~ x because 
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I A — n 2 \ > n. If n < k < 2n, then |A — k 2 \ > k 2 — n 2 — n > n, so 



k a /\X — k 2 \ < 2 a n a 1 ; finally, if k > 2n then n < k/2, and therefore, 

(2.12) |A - k 2 \ > k 2 - n 2 - n > k 2 - (k/2) 2 - k/2 > k 2 /2, 

so k a /\X - k 2 \ < 2k a ~ 2 < 2n a ~ 2 because a< 2. Hence (EH) holds. 

Next we consider the case where n < \t\ < n 2 . Since |A — k 2 \ > \t\ 
we have, for k < 2n, that k a /\X - k 2 \ < (2n) a /\t\. If k > 2n then we 
obtain, as above, that (12.121) holds, thus 



k a /\X - k 2 \ < 2k a /k 2 < 2n a /n 2 < 2n a /\t\, 

which proves (|2.9|) . 

Consider now the case where \t\ > n 2 . If k 2 < 4|t| then ( since 
|A- k 2 \ > \t\) 

k a /\X-k 2 \ < k a /\t\ < 2 a \t\ a/2 /\t\. 

If k 2 > 4\t\ > An 2 then (jTT^ holds, thus k a /\X - k 2 \ < 2k a - 2 < 
2|4t|( a - 2 )/ 2 ; which completes the proof of Lemma El 

□ 

3. By Proposition 0J for each k there is a disk of radius Rk{o) = 
k x ~ a / (8M) with the property that the operator L + zB has exactly one 
simple eigenvalue E k (z) in the strip H k . If a G [0, 1), then R k (a) f oo 
as k — >• co. 

Let us fix a 6 [0, 1) and n G N. If m > n then A n C A m , so for each 
z G A n 



Sp(L + zB)n ( |J F m j C |J tf n 

\m>n / m>n 



where K m is defined in (|2.7j) . Set 

W n = {A G C : -n < i?e A < n 2 + n, | Jm A| < n}. 

Proposition 6. Under the conditions \2.1)) . and K2.3j) . if a G 

[0, 1), i/ien /or eac/i z G A n 

(2.13) Sp(L + zB) C W n U [J AT m . 
Moreover, the projector 

(2.14) P,(z) = — [ {X-L-zB)- 1 dX 

2m Jdw n 

is well defined for z G A n , and 

(2.15) dimP*(z) = n. 
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Proof. Set H = {A G C : Re\<n 2 + n}. Then 

(2.16) sup k —. = n a ~ l for A G H \ W n . 

A — k 2 



k 



Indeed: if k < n, then |A — /c 2 | > n, so fc a /|A — A; 2 | < n a 1 ; if A; > n, 
then |A - k 2 \ > k, thus fc a /|A - k 2 \ < k a ~ l < n 01 " 1 because a G [0, 1). 
By (BUSD and (f2~TTTl. we obtain that if \z\ < n l - a /&M, then 

\z\ ■ \\BRl\\ < 1/2 for X e H \ W n . 

Therefore, in view of ()2.6|) . for each z G A n , 

+ zS) n (H \ W n ) = 0, 

which proves ()2.1Hj) because C = H U U m >n 
Moreover, the projector 

P*0) = / (A - L - zB)- l d\. 
2m JdW n 

is well defined for each z G A n , and since its dimension is a constant, 
we obtain that dim P*(z) = dimP*(0) = n. □ 

3. The Taylor coefficients of analytic functions E n (z). 
1. For each n G N, we consider the rectangles 

(3.1) II = n(n,s) = {A G C : |i?e(A-n 2 )| < n, |JmA| < s}. 
Then the one-dimensional Riesz projector 

(3.2) Pn ( z ) = -L [ (\-L-zB)- l d\ 

is well defined for \z\ < R n and does not depend on s for s > n + 1 as 
it follows from ()2.13j) and (J2.7)) . The integrand in ()3.2)1 is an analytic 
function of (A, z) G (if n \ II) x A n . 
Since 

(3.3) E n (z)P n (z) = ^- [ \{\-L-zBY l d\ 

2?u Jan 

we obtain that 

(3.4) EJz) = Trace (— [ A(A — L — zB^dX 

V 27 ™ J an 

The formulas ()3.2j) - ()3.3|) are basic for what follows in this section. 
They are used to derive formulas for the Taylor coefficients of E n (z), 
and to obtain a trace formula. 
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Let 

oo 

(3.5) E n (z) = S ^2a k (n)z k 1 a {n) = n 2 , 

k=0 

be the Taylor expansion of E n (z) at 0. 

Proposition 7. Under the conditions \2.1\) and \2.2\l with a G [0, 
we have: 

(3.6) a k (n) = £ JL \(R x (BR° x ) k e 3 , ej )dX, 
where 

[ \(R° x (BR° x )%,e j )dX = if \j-n\>k ] 
Jan 

(3.7) a k {n)= £ ±-l (A - n 2 ){R\{BR\fe h e )d\- 

\j-n\<k m JdIi 

(3.8) a k (n) = for odd k; 



(3.9) |a fc (n)|<2(2fc + l) J 1 4 _^ 1 , k>2. 
Proof. By Q and (P, 

- « OO OO 

(3.10) P B (z) = — / £i$(BJ$)***dA = E^ n ^ 

711 ^ 9n k=0 k=0 

where the integrand-series converges absolutely and uniformly for z 
A n and A G 911, and 

(3.11) Pfc ( n ) = J_/" R° x (BR° x ) k d\, fc = 0,l,2,... 

27TZ y OT 

are the Taylor coefficients of P n (z) G (|3.2|) . 
We have 

(3.12) Po(n)e„ = e n , Po(n)ej = for j 7^ n. 
Moreover, for each = 1, 2, . . . , 

(3.13) Pk(n)ej = if |j — n\ > k. 
Indeed, 



(3.14) Pk {n)e 3 = -L /" R^BR^dX. 

2 ™ J an 
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Since Be u is a linear combination of e v -\ and e u +\, while R\e v = 
e„, the singularity (or its power) could appear in the 



A — v 2 A — n 2 

integrand only if \j — n\ < k. Therefore, if \j — n\ > k, then the inte 

grand is an analytic function on IT, so the integral vanishes. 

Since dim P n {z) = 1, 

(3.15) ^(P B (z) Ci , ei > = l, 

3 

which implies, in view of (|3.12j) and ()3.13|) . that 

(3.16) ^(po^e^ej) = 1, 



(3.17) ^(p fe (n)e i)e ,) = 0, k = 1,2, — 

i 

Set 

oo 

(3.18) E n (z)P n (z) = Y,dk(n)z k . 

k=0 

Then, by (Q and (EHUl) . 

(3.19) d fe (n) = y^fl y (n)pt-i;(B). 

i/=0 

Now dSHH) and (j3~T7j) imply, in view of HH) and dSH), that 

(3.20) Ofe(n) = 2J (4(n)e 3 -,ej). 

By f!3.3|) . taking into account the power series expansion of the resol- 
vent, we obtain 

(3.21) a k (n)= £ -L / X(R° x (BR° x ) k e j , e^dX 

\j-n\<k m Jm 

Since is a diagonal operator, and -Be., is a linear combination of e^-i 
and ej + i, we have 

(3.22) X{(BR° x ) k ej , e 5 ) = Vj if k is odd. 

Therefore, the same argument that explains (|3.13|) (see (J3.14|) and the 
text after it) shows that the integrals in ()3.21|) are equal to zero if 
\j — n\ > k, which proves ()3.6|) . 

Since Bej is a linear combination of and e^+i and R x is a di- 
agonal operator, we obtain for odd k that (BR x ) k ej is a finite linear 
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combination of vectors e„ such that v — j is odd number, so v ^ j. 
Therefore, if k is odd, then for each j the integrands in (J3.6)) are equal 
to zero, which proves (|3.8|) . 

By dSHSD , (jSH and (|3~T7j) we have 

|j-n|<fc 

thus (ESI implies (pTTjl . 

Next we prove ()3.9|) . Let us replace the contour 911 in ()3.7|) by the 
circle (9i^ n = {A : |A — n 2 \ = n}. Fix j with |j — n\ < k and consider 
the corresponding integral. The integrand does not exceed 

sup (|A-n 2 | • \\R° X \\ ■ \\BR° x \\ k ) . 

\GdK n 

By (|23J), we have, for a G [0, 1), 

\\BR° X \\ < AMn*- 1 if A G dK n . 
On the other hand |A — n 2 \ = n on 9iif n , and 

llirfll = sup — — < - for AG dK n . 
j X-f n 

Thus, for each j, the integrand does not exceed 2(4M) fc n( Q_1 ' )fc and the 
length of dK n is equal to 2-nn, which leads to the estimate ()3.9|) . □ 

2. Next we give another integral representation of the coefficients 
a k (n). 

Proposition 8. Under the conditions \2.1\) and \2.2\) with a G [0, 1) 
we have, for each k > 2, 

(3.23) a k (l) = y?fc(l), a k (ri) = (p k (n) - ip k (n - 1), n > 2, 

(3.24) ^(n) = E i / A(^(5^) fc e,, e^A, 

(3.25) /i„ = {AGC: i?eA = n 2 + n}, 
and 



/ X(R° x (BR° x ) k e J ,e J }dX = if \j 



n\ > k. 



Moreover, 

(3.26) \Mn)\ < ^ {1 %_ 2 , k>l, C k = (2k + l)(8M) k . 



n 
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Proof. Letting s — > oo in (jH.fi j) we obtain (j3.23j) - (j3.25j) . To justify this 
limit procedure, we have to explain that 

(i) the integrals over h n and /i n _i converge; 

(ii) the integrals over horizontal sides of dH(n, s) G (13. lj) go to zero 
as n — > oo; 

(iii) the integrals over h n are equal to zero if \j — n\ > k; 

(iv) ajfc(l) = ^fc(l)- 

Indeed, (i) and (ii) hold because the integrand in (j3.6j) . for each even 
k > 2, is a linear combination of rational functions of the form 

(3 ' 27) QiJ ' x> = (>-m-ffl--o-Ay J=u °- 

and therefore, the integrand decays faster than 1/|A| 2 as |A| — > oo. 

(iii) If j — n > k (respectively n — j > k), then the integrand is a 
sum of terms (j3.27j) with j , . . . , j k > n (respectively jo, ■ ■ ■ ,jk < n )• 
Consider the contour that consist of the segment {A G h n : \Im A| < s} 
and the left half (respectively right half) of the circle with center n 2 + n 
and radius s. Since the integrand is an analytic function inside the 
contour, the integral is equal to zero. Letting s — > oo we obtain that 
the integral over h n is zero, because the integral over the half-circle 
goes to zero due to the fact that the integrand decays as 1/|A| 2 or more 
rapidly. 

The same argument shows, for each j, that the integral over the 
imaginary line ReX = equals zero, which explains (iv). 

Finally, we prove (j3.26j) . By (j3.24j) . the function ipk{n) is a sum of 
at most 2k + 1 integrals over h n of the form 

(3.28) -L / X(R x (BR° x ) k e 3 ,e 3 )dX. 
2m Jh n 

The absolute value of the integral (j3.28j) does not exceed 

(3.29) — / F(t)dt, where F(t) = \\XR x (BR° x ) k \\, X = n 2 +n+it. 
2tt J r 

Next we estimate from above F(t) < \\XR X \\ ■ \\BR x \\ k . Lemma El gives 
estimates of the norm on each of the three sets 

h = {t: \t\ < n}, h = {t: n< \t\ < n 2 }, h = {t : \t\ > n 2 }. 

On the other hand we have 



n + 1, t G h 

(3.30) ||A^|| = ' < { 2n 2 /\t\, tel 2 



\n 2 + n + it\ 



\n + it\ 



tel. 
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If we combine ()3.3()j) with the estimates ()2.8|) - (j2.1()|) from Lemma Owe 
get 

2n (4Mn a - 1 ) k , t G h 
(3.31) F{t)<{ 2n 2 (AMn a \t\- l ) k , t G J 2 



2(8M|t|( a " 2 )/ 2 ) fc , te/ 3 



Therefore, since 



/ = / F(t) + I F(t) + [ F(t), 

JR J h Jh Jh 



the estimates ()3.3H) imply that ()3.26|) holds. 



□ 



3. The formulas ()3.23j) and ()3.24j) could be used to find the Taylor 
coefficients of E n (z). Indeed, under the conditions (|2.1|) and ()2.2|) with 
a G [0, 1), a computation based on the standard residue approach shows 
that 



(3.32) cp 2 (n) 



2n + l' 



/o o<->\ / \ b n c n b n b n ^\C n c n -^-i b n b n —iC n c n - 

(3.33) fiyn) — 



(2n+l) 3 (2n + l) 2 (4n + 4) 4n(2n + l) 2 ' 

For any off-diagonal sequences b,c G (|2.1|) + (|2.2j) . it follows from 
()3.26|) that as n — > oo 

(3.34) <p k (n) ^ if a < 2/3, fc > 6, 
and by (jS32), (HO^ 

(3.35) ip 2 {n) -»• 0, if a < 1/2, y> 4 (n) -»• 0, if a < 3/4. 
Now, by (glMD and (JOS), 

oo 

(3.36) ^ afc (n) = if k > 6, a G [0,1/2]. 

n=l 

If (OJ) holds, then 

(3.37) <p 2 (n) = 



n 2a 



(3.38) p 4 (n) 



2n + V 



(2n + l) 3 (2n + l) 2 (4n + 4) 4n(2n + l) 2 
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By Q3I23D and (f3~3Tj) we obtain 
(3.39) 

a2 (i) = V2 (i ) = --, a2W = L_L___ for „> 2 . 

Observe that y> 2 (n) -> if a G [0, 1/2), while y» 2 (n) -> -1/2 if a = 1/2. 
Thus we have 

oo 



(3.40) = for a G [0, 1/2) 

71=1 

and 

(3.41) ^a 2 ( n ) = -- if a = 1/2. 

n=l 

By fHOH]) we obtain that ^4(71) — > if a G [0, 1/2], so (jSHSD yields 

CO 

(3.42) 5^«4H = if a G [0, 1/2]. 

71=1 

4. ASYMPTOTICS OF E n (z) 

In this section we study the asymptotic behavior of E n (z) for large n. 
Our approach is based on the fact that the eigenvalue function E n (z) 
satisfies a quasi-linear equation. Of course, the same estimates and 
formulas could be found if one follows the Raleigh-Schrodinger scheme 
with recurrences for the Taylor coefficients 

oo 

K z ) = ^2a2k(n)z 2k , a (n) = n 2 , 

k=0 
oo 

/(z)=X>2>, fj G ^(N), f = e n , 

as they would come if one substitute the above formulas into (|4.1j) . 

1. Throughout this section we assume that (|2.1|) and (|2.2j) with 
a G [0, 1/2] hold, but after (|4.1U|) we assume that (|1.3|) holds also. 

Suppose that n and 2; G A n are fixed and A = E n (z) is the corre- 
sponding eigenvalue of the operator L + zB. Then we have 

(4.1) (L + zB)f = \f 

for some / ^ 0. Let P be the projector defined by Px = (x,e n )e n , and 
let Q = 1 — P. The equation (|4.1|) is equivalent to the system of two 
equations 

(4.2) {\-L)f 1 = zPB{f l + f 2 ) 1 



14 



PLAMEN DJAKOV AND BORIS MITYAGIN 



(4.3) (\-L)f 2 = zQB(fi+f 2 ), 

where fx = Pf, f 2 = Qf ■ The operator A — L is invertible on the range 
of the projector Q; we set 

(4.4) De k = yz^Zk if k ^ n, De n = 0. 

Then D is well defined in £ 2 , and (A — L)Dx = x on the range of Q. 
Acting on both sides of ()4.3|) by the operator BD we obtain 

(4.5) Bf 2 = zTBfi + zTBf 2 , 
where 

(4.6) T = BD. 

The operator 1 — zT is invertible for each z G A n . Indeed, since Tet = 
BR\ek foikj^n and Te n = 0, the proof of (|2.8|) shows that 

(4.7) ||T|| < AM ■ n"" 1 for A G H n . 
Thus we have 

\\zT\\ < \z\ ■ \\T\\ < 1 

for each z G A n and each A G H n . 
Solving ()4.5|) for Bf 2 we obtain 

(4.8) Bf 2 = z(l-zT)- 1 TBf 1 . 
Inserted into ()4.2|) . this leads to 

(A - L)h = zPBh + z 2 P(l - zT)~ x TBf\, 
which implies (since 1 + zT(l — zT)^ 1 = (1 — zT)^ 1 ) 

(4.9) (A - L)h = zP{\ - zT)~ 1 Bf 1 , 

where fi = const ■ e n ^ (otherwise, by (|4.3J) it follows that f 2 = 0, 
so / = fx + f 2 = 0, which contradicts / ^ 0). Since Le n = n 2 e n , the 
equation ()4.9|) is equivalent to 

(4.10) \-n 2 = z({l-zT)- 1 Be n ,e n ). 
Since 

(4.11) Be k = (k-l) a e k ^ + k a e k+1 , 
we have, by (jOll and (IP1) . that 



(4.12) Te fc = ((A; - l) a e k _x + fc Q e fc+1 ) , Te n = 0, 
and therefore, 

(4.13) (T 2k Be n , e n ) = 0, fc = 0,1,2,.... 
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Let 



Em 



n 2 + ai(n)z + a2(n)z 2 



be the Taylor expansion of E n (z). Set for convenience 

(4.14) C(z) = E n (z) - n 2 = a x [n)z + a 2 (n)z 2 + ■■■ . 
Then, by 

(4.15) ((z) = (TBe n , e n )z 2 + (T 3 Be n , e n )z A + (T 5 Be n , e n )z 6 + ■■■ 
where, by ()4.12jl . the operator T depends rationally on A = E n (z) 



C(*) 



n 



] . It is easy to see, by induction, that (|4.15|) yields a 2 k+V 



n 



0, k G N (in fact, we know this from Section 3.1, see (|3.8jl ). Thus we 
have 

(4.16) ((z) =a 2 (n)z 2 + a 4 (n)z 4 + ■■■ . 

One may use ()4.15|) to compute the Taylor coefficients of C(^)- m_ 
deed, the right side of ()4. 15|) is a power series in z which coefficients are 
rational functions of A = (+n 2 without a singularity at 0. So, replacing 
these rational functions with their power series expansions at 0, and 
replacing ( with its power expansion (J4.14)) . we obtain (comparing the 
resulting power series expansion on the left and on the right) a system 
of equations for the coefficients a 2 (n), 0,4(71), 

Next we compute some of these coefficients. By (J4.11|) - (J4.15J) it 
follows that 



(4.17) 



C 



z 



1 



,2a 



n 



2d 



(n 



2n- 
l) 2a (n-2) 2a 



C 



2n + 1 -C, 
n 2a (n + l) 2a 



(2n - 1 + C) 2 (4n - 4 + C) 

(n ~ 1) 2Q i C 
2n - 1 V 2n-l 

(n - l) 2a {n - 2) 2a 



(2n + 1 -()2(4n + 4 + C) 

2a 



+ 



n 



2n + 1 



c 



(2n - l) 2 (4n - 4) 
n 2Q (n + 1) 2q 



2n 



2n + 1 

C 

An -4 



+ 



+ 



+ 



(2n + l) 2 (4n + 4) 
Hence we obtain 



1 



c 



2n + 1 



1 + 



4ra + 4 



+ ■ 



+ 



(4.18) 



02(0?, n) 



n 



\2a 



11 



2d 



2n 



2n + 1' 



n > 2; 



(4.19) 



04(0, n) 



-a 2 (n)) 



\2o 



+ 



11 



2<> 



(2n- l) 2 (2n + l) 2 



+ 
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(n - l) 2a {n - 2) 2a n 2a (n+l) 2a 



n>3. 



(2n - l) 2 (4n - 4) (2n + l) 2 (4n + 4) ' 
The same method gives 

(4.20) a 6 (ct, n) = a\(n) — a 2 (n)cr 2 (n) — a 4 (a, n)o~ 3 (n), n > 4, 
where 

(n-l) 2 >-2) 4Q n 2 > + l) 4 " 

{ ' ai{n) (2n - l) 3 (4n - 4) 2 (2n + l) 3 (4n + 4) 2 

(n - l) 2a (n - 2) 2a (n - 3) 2q n 2a {n + l) 2o (n + 2) 2a 



(2n - l) 2 (4n - 4) 2 (6n - 9) (2n + l) 2 (4n + 4) 2 (6n + 9) ' 

/ \ (n-l) 2a (n-2) 2a (2 1 \ 

(422) * 2 W ^ (2,- l) 2 (4n- 4) (2^1 + 4^4 j + 

n 2a (n + l) 2Q / 2 + 1 \ n 2a (n - l) 2a 



(2n + l) 2 (4n + 4) \2n + 1 4n + 4 J (2n + l) 3 (2n - 1) 
(4 . 23 ) *(») = &^£ + " 2 ° 



(2n-l) 2 (2n + l) 2 ' 

Of course, the case of small n requires a special treatment. For 
example, if n — 1, then with 

C = a 2 (l)^ 2 + a 4 (l)z 4 + a 6 (l)^ 6 + • • • 

we have 



c _ 3 ; ^ (c _ 3 )2 (c _ 8) 

(C - 3) 3 (C - 8) 2 + (C - 3) 2 (C - 8) 2 (C - 15), ' + ' 
which leads to 

1 1 2 2a 

(4.24) a 2 (l) = --, a 4 (l) = - - — 

(compare with (ET%ty . lCT37|) . and 



(425) ae(1) " "3^ ~ 3^8^ + 3^ " 3*" 

2. The following lemma gives the asymptotic behavior of a 2 fc(a,n) 
as n — > 00. 

Lemma 9. Under the condition M.ty) . if a & [0, 1), i/ien 
(4.26) a 2fc (a,n) = 0(n 2fe ( a - 1 )). 
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Proof. We prove (|4.26|) by induction in k. If k = 1, then (j4.18j) yields 
(4.27) a 2 (a,n) = 0(n 2 ^). 

If k = 2, then ()4.19|) gives a^a, n) as a sum of two expressions. For 
the first one we obtain, in view of ()4.27|) . that 

*<«•"> (H^S + (STi?) = (« !, "- ,, »'°(« 2( °- I, » - °<" 4, "- I, >- 

The remaining part of ()4.19|) is 

(n - l) 2Q (n - 2f a n 2a {n + l) 2a 



^ 4 ' 28 ^ (2n- l) 2 (4n-4) (2n + l) 2 (4n + 4) ' 

Each term of this difference is 0(n Aa ~ 3 ). But (jj^HD is 0(n^ a ^) due to 
the Mean Value Theorem. Indeed, let f(m) = m 2a (m + l) 2a (4m + 4) _1 
and g(m) = (2m + 1)~ 2 . Then f!4.28|) may be written as 

f(n-2)g(n-l)-f(n)g(n) = (f( n -2)-f(n))g(n-l)+f(n)(g(n-l)-g(n)). 

Since 

f(t) = 0(f(n)/n), g'(t) = 0{g{n)/n), for te[n-2,n], 

by the Mean Value Theorem the expression (J4.28JI is 0(n 4a ~ A ) which 
proves (|4.26J) for k = 2. 

Fix k > 3 and assume that (|4.26|) holds for 1, . . . , k — 1. Then by 
(ETToT) and (jUm we obtain, in view of ()4.17jh that 



o-2k — (T 2k 1 Be n ,e n ) + ^C row „, //o" • • • </ 



m fe _i u '2 ' ' ' U 2(fc-1)' 



where mi + 2m 2 + ••• + (& — l)mfc_i = fc, and 

1 

u," — v 

In addition, for each term of the sum, we have 



T = BD with De n = 0, De v = — -e v . 



a n m i... n m *-i _ n(n 2k ( a - 1 h 

L'mi...m t _ 1 "2 U 2(fc-1) ~~ W \ U )• 

(See ()4.2U|) — f)4.23f) for the case k — 3.) Thus Lemma El will be proved if 
we show that 

(4.29) (T 2k - x Be n , e n ) = 0(n 2k ^). 
Set 

(4.30) B = B + i + B_i and T = T +1 + T_ 1; 
where 

(4.31) 5 +1 e fc = fc a e fc+1 , £Ue fc = (fc - l) a e fc _ 1; 
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and 

(4.32) T +1 = B +1 D, r_! = B^D. 
Then 

(4.33) (T 2k - 1 Be n ,e n ) = Y,^)i 

where the summation is over all 2fc-tuples e — (ei, . . . , E2k) with e v = 
±1, and 

(4-34) u(e) = (T £2k _ 1 ---T £2 B £l e n ,e n ). 

Let 

(4.35) 5(e) = (5i, 5 2k ), S v = 5 u (e) = eH h£„, i/ = 1, . . . , 2k; 

then T £v ■ ■ -T £2 B ei e n = const- e n+ $ u . Therefore, since De n = 0, we have 
u(e) 7^ if and only if 5 2 k = and 5 U ^ for v ^ 2k. 
Now iffTSSjl implies that 

(4.36) (T 2k ~ l Be ni e n ) = £ + u;(-e)], 

where the summation is over the set e + of all 2/c-tuples e such that 
S n u(e) > for z/ = 1, . . . , 2k — 1. Since the cardinality of e + does not 
exceed 2 2k , (J4.29)) will be proved if we show, for each e G e + , that 

(4.37) w(e) + w(-e) = 0(n 2fc(a - 1} ). 
By dOH-dOl we obtain 

lL=i (^(2n + ^)) 

Now, as above, the Mean Value Theorem may be used to show that 
()4.37|) holds. This completes the proof of Lemma M 

□ 

3. Proof of Theorem^ By Proposition |U we know that, with R n = 
n 1-o /(8A0> 

|£„(z)-ra 2 |<ri for zeA n = {(: \(\<R n }. 
Lemma 10. For each k = 1, 2, . . . , 

(4.38) \a k (n)\ = ^ |^ fc) (0)| < (8M) V- fc(1 - a) . 

k\ 1 1 

Proof. Indeed, E n (z) is analytic in A n . Therefore, the Cauchy inequal- 
ity for the Taylor coefficients of E n (z) at gives f)4.38j) . □ 
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Now, for \z\ < R, we obtain 



n 13-14a ' 



(4.39) \E n (z) - n 2 - ]>> 2 ^ 2fc | < £ \a 2k (n)\R 2k < 

k=l k=7 

where C n = (8MR) U £ fc > (8Mi?) 2fc /n 2 ^ 1 "^ is a bounded sequence. 

On the other hand (|4"Tg|) and (}4~T9"|) imply that 
(4.40) 

. , (l-2a) (a 2 -a) (1 - 2a)(8a 2 - 14a + 3) 2a _ 5 . 
a 2 (a, n) = - — — - + , n + - - — ^ +0(n 2a 5 ) 

and 

(4.41) a 4 (a,n) = 0(n 4a " 6 ). 
The formulas (j33nj) -(|4~23 |) yield 

(4.42) a 6 (a,n) = 0{n 6a ' w ). 
Analogous computations show that 

(4.43) a 8 (a,n) = 0(n 8a - 14 ). 
Finally, by Lemma El 

(4.44) a 10 (a, n) = O(n 10a - 10 ), a 12 (a, n) = 0(n 12a - 12 ). 

Now (jCTjl - (j04]l imply (fTKjl . Indeed, if a e [0, 1/2], then 2a - 5 > 
4a — 6; moreover, 

12a - 12 < 10a - 10 < 2a - 5 

and 14a - 13 < 2a - 5, thus (jEEJ) holds. 

If a G [1/2, 2/3], then 2a — 5 < 4a — 6; so, since 

12a - 12 < 10a - 10 < 4a - 6 

and 14a — 13 < 4a — 6, we obtain that (jl.8|) holds. This completes the 
proof of Theorem El 

4. We consider separately the case where a = 1/2 in the following 
theorem. 

Theorem 11. If \z\ < R, then 

2 n 2 , o 4 

(4.45) E n (l/2, Z )=n 2 - — - + {l/n% 
Proof. Ha = 1/2, then (EL"3HJ) implies 

6 OO 

(4.46) \E n (z) - n 2 -^a 2k z 2k \ < £ |a 2fc (n)|ii! 2fe < -J, 

fc=l k=7 
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where C n = (8Mi?) 14 ^ k>0 {9>MR) 2k /n k is a bounded sequence. On 
the other hand, from (gj^ - jg^ it follows that 

(4.47) a 2 ' 



(4.48) a 4 



4n 2 - 1 ' 
1 1 



4(2n + l) 3 4(2n-l) 



3 



(4.49) a 6 = --— -— —i— — -4 ' 



(2n + 3)(2n + l) 5 (2n - 1) (2n + l)(2n - l) 5 (2n - 3) ' 

The same approach that leads to ()4.18|) - (j4.23|) gives 

-327- 16080n 2 - 63136n 4 + 29440n 6 + 39168n 8 

^ ' ' a& ~ 32(n - l)(n + l)(2n - 3)(2n + 3)(2n - l) 7 (2n + l) 7 ' 

and 
(4.51) 

3915 + 280676n 2 + 2496992n 4 + 2635904n 6 - 3111168n 8 - 1158144n 10 
° 10 ~ 8(ra - l)(ra + l)(2n - 3)(2n + 3)(2n - 5)(2ra + 5)(2n - l) 9 (2n + l) 9 
By (IQ81 - (I43T|) we obtain 

(4.52) a2( i/2,n) = -^-^L + 0(n- 6 ), 

(4.53) a4 (i/ 2 ,n) = -^ + 0(n- 6 ), 

and 
(4.54) 

a 6 (l/2, n) = 0(n~ 8 ), a 8 (l/2, n) = 0(n- 10 ), a 10 (l/2, n) = 0(rT 14 ). 
In addition, Lemma H3 implies that 

(4.55) a 12 (l/2,n) = 0(n~ 6 ). 
Now (IQoT) follows from (l4~47)|) and (14321 - d4~H5l . 

□ 

Remark 12. 

We evaluate ai2 in ()4.55j) by using the general estimate (|4.2fij) from 
Lemma HU However, a direct computation of the coefficients a2k('n) 
for 6 < k < 14 shows that each of them is 0(l/n 16 ). Estimating the 
remainder as in the proof of Theorem ^2 we get 

(4.56) a ik(n)z k = 0(l/n 14 ), \z\ < R. 

k>15 
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So, by (|4.51jl . we have 
(4.57) 

E n (z) = n 2 + a 2 z 2 + a 4 {n)z /L + a 6 {n)z 6 + a 8 {n)z 8 + 0{l/n u ), \z\ < R. 
It follows from here, in view of ()4.47JI - ()4.5U|) . that 

6 1 

(4.58) E n (z) =n 2 + J2 + 0(l/n 14 ), 



rr 

k=l 



where 

(4.59) P 1 ( Z ) = -- P 2 (z) = P s (z) 



4 ' 32 JV ' 64 



, , -2z 2 - 21^ 4 + 28z 6 „ , , -8z 2 - 144^ 4 + 1920z 6 + 153^ 8 
v ; 512 ' 5V ; 8192 

„ . x -2z 2 - 55z A + 5192z 6 + 880z 8 

bV ; 8192 
See further discussion in Section 7.3. 

5. Analytic continuation of eigenvalues and Regularized 

Trace 

1. Each eigenvalue E^z), as we have seen in Proposition HJ is well 
defined and simple if \z\ < Rk = k l ~ a /8M. We are going to show 
that it is possible to continue Ek(z) analytically as z is moving along 
a smooth curve which goes around singular points ( G S, where S is a 
countable set without a finite point of accumulation. 

Fix n6N and consider the rectangle 

W = W n = {X E C : -n< ReX<n 2 + n, \ Im X\ < n}. 
By Proposition El the projector 



P*(z) = — [ (A - L - zB)- 1 d\ 

27fi Jg w 



is well defined for z 6 A n and 

dim P*(z) = n. 
Consider the analytic functions 

(5.1) (Tj(z) = Trace ( — [ X J (A - L - zB)~ x dX 1 , 1 < j < n. 



2m 



dW 



If \z\ is small, say \z\ < e < Ri, then 

n 

(5.2) a J (z) = J2(Ek(z)y, l<j<n, 



k=l 
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where all Ek(z) are well defined. Moreover, 

n n 

(5.3) ]J(\-E k ) = J2Qn- j (E)\ j , 

1 

where {Qi}™, Qo = 1, are symmetric polynomials of {E k }'1. But {<Jj}i 
is a basis system of symmetric polynomials (see, e.g. [17]). and there- 
fore, 

(5.4) Qj = q 3 (a) 
are polynomials of cr's. Thus 

n n 

(5.5) fl(A - £,0 = J>-i(*C*))A', 

1 o 

at least for small \z\, say \z\ < e. However, the coefficients Cj(z) = 
q n -j(a(z)) are well-defined by ([5.4]) . ([5.1]) in the entire disk A n and 
analytic there. The factorization ([5.3]) becomes 

n n 

(5.6) n(A-^)) = 53c,(^ 1 

1 o 

and the equation 

n 

(5.7) c(z, A) := Cj{z)X j = 0, |z| < i? n , 

o 

defines over A n the surface 

(5.8) G n = {(A, z) eC x A n : c(z,\)=0} 

with n sheets and possible branching points if the polynomial Ylo Cj(z*)X 
has multiple roots. Such a point root of the resultant 

(5.9) r(z) = R(c(z,-),c' x (z,-)) 

of the polynomial c(z,X) and its derivative c' A . Notice that r(z) is an 
analytic function of z, \z\ < R n , because the resultant is a polynomial 
of Cj(z) eflEU). If z = 0, then 

(5.10) ^(OjA^nCA-* 2 ), 

o 1 

and all zeros are simple. Therefore, r(0) 7^ 0, so the resultant r(z) is 
not identically zero. Thus the set 

(5.11) S n = {z e A n : r(z) = 0} 
is finite. By Proposition 0] we can conclude that 

(5.12) S n C S n+ i and S n+i n A n = E n . 



TRACE FORMULA 



23 



Thus the set 

(5.13) S = |j£„ 

is countable and has no finite points of accumulation. 
We have proved the following. 

Proposition 13. Under the conditions of Proposition ^| there is a 
countable set S without finite accumulation points such that if 

7 = {z(t) : < t < T}, z{0) = 0, 7 nS = 

is a smooth curve then each eigenvalue function E^z), Ek(0) = k 2 , 
can be extended analytically along the curve 7. 

2. We define Spectral Riemann Surface (SRS) of the pair (L, B) as 

(5.14) G = {(X,z)eC 2 : (L + zB)f = \f, f G e (N), / ^ 0}. 

Proposition 14. Under the conditions of Proposition^ for each z G" S 
the surface G has infinitely many sheets over a neighborhood U £ 3 z for 
small enough e(z) > 0. Each branching point z* e S is of finite order. 

Proof. Everything has been already explained. The surface G over A n 
is defined by f|5.7|) . and, by ()5.7|) - ()5.11|) . X(z) has branching points 
z* G A n of order < n. □ 

3. We follow the 1975 Schafke construction (see [20], pp. 88-89), 
as it is presented by H. Volkmer [10], to analyze whether the Spectral 
Riemann Surface G is irreducible. 

Let k,j e N. We call k and j equivalent, k ~ j, if there is a smooth 
curve 

ip: [0,T]-C\£, <p{0) = <p{T) = 0, 
such that the analytic continuation of Ej.{z) along ip leads to Ej(z). (A 
Spectral Riemann Surface G is irreducible if N is the only equivalence 
class, i.e., k ~ j for any k,j e N.) 

Such construction, carried for each fceN, defines a mapping 

such that TT v -i(j) = k, where ip _1 (t) = y(T — t). With R n — > 00, we 
have for some n that max[ 0i Ti 1^(^)1 ^ ^n- Therefore, by Proposition 

(5.15) n<p(k) = k if k > n. 

Lemma 15. Let M. be an equivalence class (or union of equivalence 
classes), and n G N. Then the function 

(5.16) E n (z)= £ 



24 



PLAMEN DJAKOV AND BORIS MITYAGIN 



(which is well-defined and analytic for small enough \z\) can be ex- 



tended analytically on the disk {z : \z\ < R n }, R r , 



n 



l-a 



/(8M). 



Proof. Take any smooth curve ip : [0, T] — > A n \ S, such that <^(0) = 
(fi{T) = 0. Then 7r = 7T V : M. — > M. is a bijection, and (|5.15|) holds, so 
7r permutes the finite set {k e M : k < n}. Therefore, E n (z) can be 
continued analytically, term by term in ()5.16|) . and the result will be 

E 

j&M,j<n 



E 



Ej(z) 



E n (z), zeA n \S, 



keM,k<n 



i.e., the same function. By Proposition EJ if \z\ < R n , then we have 
exactly n eigenvalues on the left of the line h n = {Re A = n 2 + n}, and 
all of them lie in the rectangle W n . Therefore, 



(5.17) 



E n {z) <n(n 2 + 2n). 



So the function E n (z) is analytic and bounded on A n \ S, while the set 
A n n S is finite. Thus, it is analytic in the disk A n . □ 

The inequality ()5.17j) cannot be improved essentially because 

n n 

^£ fc (0) = ^2k 2 = n(n + l)(2n + l)/6 ~ n 3 . 



However, we can regularize Ek(z) by considering Ek(z) — Ek(0), where 
E k (0) is real. 

Again by Proposition |U1 if \z\ < R n , then the operator L + zB has 
n eigenvalues that lie in the rectangle W n , so the absolute value of the 
imaginary part of each of these eigenvalues is less than n. Therefore, 



(5.1? 



Im[E n (z)-E n (0) 



< n 2 . 



By Borel-Caratheodory theorem (see Titchmarsh [30], Ch.5, 5.5 and 
5.51), if g(z) is analytic in the disk \z\ < R, g(0) = and \Img(z)\ < C, 
then \g(z)\ < 2C for \z\ < R/2. Thus (|OH|l implies 



(5.19) 



E n {z) - E n (0) < 2n 2 , for \z\ < R n /2. 



This conclusion is valid for each equivalence class, or union of equiv- 
alence classes Ai; in particular, for M. = N. 

4. Definition of the regularized trace tr(z). Now we are ready to 
define an entire function tr(z), the regularized trace of L + zB, under 
the conditions (EHJ and (Q with a < 1/2, or flO)) with a = 1/2. 
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For small z, \z\ < R\ = 1/(8M), all E n (z) are well defined, and 



oo / oo 

2fc 



(5.20) tr(z) = ( E ni z ) ~ O = E I E a ^( n ) z 

n=l n=l \fc=l 

00/ 00 

= ^2 a 2( n ) z2 + a 4 (n)z 4 + ^2 a 2k{n)z 



2 k 

n=l \ k=3 



z 2 



lim <^2(p) + -2 4 • lim <Pi(p) + 7 7 

p— »00 p— >0O ' ^ ' ^ 

n=l 3 



By ()3.32|) - (j3.35p . the latter limits are well defined and the third term 
is an absolutely convergent series. Indeed, by (j4.38j) . Lemma ITTH we 
have for \z\ < 1/(8M) that 

e(em»)i)ni^e(eJS)^ 

n=l \fc=3 / n=l \k=3 / V ' 

00 

E—ttz — r • < 00 for a < 2/3. 

n=l 

Therefore (|5.2U|) defines tr(z) as an analytic function in the disk \z\ < 
1/(8M). 

Fix N eN and consider the analytic function En(z), z G An, given 
by Lemma [TBI in the case where Ai = N. For small \z\ we have 

00 

(5.21) tr{z) = E N (z) - E N {0) + {E n {z)-n 2 ). 

n=N+l 

The same formula gives the analytic extension of tr(z) on A^v because 
the series on the right side of (J5.21|) converges uniformly on A^v- Indeed, 
with E n (z) = a2{n)z 2 + a 4 (n)z 4 + • • • , we have 
(5.22) 

00 / 00 \ / 00 \ 0000 

(E n { z ) - n 2 ) = Ij2 a ^ n ) ) z 2 + E a ^ p 4+ E E a ^( n > 2k - 

N+l \N+1 / \N+1 J n=N+l k=3 

By JOJ) we obtain, for n > N + 1 and \z\ < R N = N^/fiM), that 
(5.23) 

J2 \a 2k (n)\ < E^tr^ 4M « ^ U 

fc=3 fc=6 ^ ' 

where C(N, a) = N^ 1 '^ Efc> 6 ( 4A; + 2 ) 2 ~ fc < 00 • Now ' in view of 

the estimate ()5.23|) implies that the series in (j5.21|) converges uniformly 
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in An if a G [0,1/2], thus tr(z) can be extended analytically in the 
disk A N . Since U^A^r = C this defines tr(z) as an entire function. 

5. Proof of Theorem^ According to the previous subsection, tr(z) is 
an entire function. Therefore, it is enough to prove (jl.7|) only for small 
\z\, or to evaluate its Taylor coefficients. By (|5.22J) tr(z) = ^2kZ 2k 
where 



Mk = VVfcO) = lim f 2 k(p)- 

1 ' p— >oo 

n=l 



If a < 1/2, then we have, by (HOI) and (l3~3o]) . that 
lim if2k{p) = A: = 1, 2, ... , 

p— »oo 

and therefore, tr(z) = 0. 

If a = 1/2, then (HOD and (EHoT) imply 

lim <f2k{p) = fc = 2, . . . , 

p^oo 

and by (|3.32j) . if the limit I = limbkC^/k exists, then 



lim (pzip) = lim 

p— *oo p— »oo 



6 C 



2p + iy 2' 

so tr(z) = — (£/2)z 2 . This completes the proof of Theorem 

6. Spectral Riemann Surfaces 

In our analysis of the regularized trace it was important to see that 
by inequality (|4.38J) from Lemma ITOl 

-fe(l-a) 



a k (n) = ±\EW(0)\<(8M) k n 1 - k < 1 



so the series 

oo 

(6.1) Y,\ E n ) (°)\< 00 if a<l-2/k. 

71=1 

Therefore, for every subset AicN, the partial sum 

(6.2) e<*>{M) = ^(0) 

is well defined. 

On the other hand, ()5.19|) and the Cauchy inequality for the Taylor 
coefficients yield 



1 



k\ 



2n 2 



(6.3) — Ei k \0) < — - = 2{16M) k n 2 - k{1 - a) 



(R n /2Y 
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So, if a < 1 — 2/k, then 



(6.4) 



lim 

n 



lim 

n 



£f(0) 

m£M,m<n 

Therefore, the following statement is true. 



E ^ fc) (o) 



0. 



Proposition 16. If a < 1 — 2/k, then we have for each equivalence 
class Ai of the Spectral Riemann Surface of the pair (L, B) that 



(6.5) 



4. Finally we show that some Spectral Riemann Surfaces are irre- 
ducible, which is the claim of Theorem El 

Proof of Theorem First we consider the case where 11. Hj) holds 
with a = 1/2. By Proposition [TBI Theorem El will be proved if we show 
that there is no proper subset M. C N with the property (|6.5j) for a 
fixed k > 2/(1 — a). Indeed, then N will be the only one equivalence 
class, which implies that the Spectral Riemann Surface is irreducible. 

If a = 1/2 then k = 6 is the least even k for which k > 2/(1 — a). 
By ()4.49j) we have 

(6.6) ^£i 6) (0) = a 6 (l/2, n) = V>(n) - if>(n - 1), n 



2,3,... 



where 
(6.7) 



i/)(n) 



1 



(2n- l)(2n+ l) 5 (2n + 3)' 
On the other hand, from (J4.24j) . with a = 1/2, it follows that 

1 



5-3 5 ' 



(6.8) a 6 (l/2,l) = V(l) = 
In view of (j6~H|) - (l6~H|l . 

(6.9) a 6 (l/2, 1) < 0, a 6 (l/2,n)>0 for n > 2, 
and 

oo 

Ea 6 (l/2,7i) = -a 6 (l/2,l). 



n=2 



Certainly, 
(6.10) 



a 6 (n) 



then .M = N. 



This proves Theorem 01 for a = 1/2. 
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If a G [0, 1/2), then 4 > 2/(1 — a), so, in view of Proposition ITHl and 
the above discussion, the Spectral Riemann Surface corresponding to 
a G [0, 1/2) will be irreducible if all but one terms of the sequence 

a 4 ( a) n) = ^(0) 

have the same sign. Below, in Lemma IT71 we show that this is true if 
a G [0,0.085] and a G [(2 - v / 2)/4, 1/2], which completes the proof of 
Theorem EJ 

5. For convenience we set 7 = 2a and 

04(7, n) = 04(7/2, n), £4(7, n) = ^4(7/2, n). 

Lemma 17. Under the above notations we have 

(6.11) a 4 ( 7 ,l) =^4(7,1) = ^"^>0, 7G[0,1]; 

(6.12) 5 4 (7,2)<0, T G[0,1]; 

(6.13) a 4 (7,n)>0 1/ 7 G [0, 0.1717], n > 3; 

(6.14) a 4 (7,n)<0 «/ 7 6 [(V2- 1)/V2,1], n > 3. 
Proof. By ()3.23j) and ()3.38j) we have that (jfi.llj) holds, and moreover, 

(6.15) a 4 (7, n) = <^ 4 (7, n) - <p±(j, n - 1), 
where 

r , ~ / n ^ 7 n> + l) 7 _ (71 - l) 7 n 7 

1 j (2n + l) 3 (2n + l) 2 (4n + 4) 4n(2n + l) 2 ' 

In particular, 

/2 27 6 7 2 7 \ / 1 2 7 
a 4 ( 7 , 2) = ^4(7, 2) - &( 7 , l)=^_-_-_j-^_-_ 

Graphing 04(7, 2) one can easily see that (J6.12)) holds. In the same way 
one can verify that the following inequalities hold: 

(6.17) 5 4 ( 7 ,m)>0 if 7 G [0, 0.1717], m = 3,4,5,6, 
and 

(6.18) fi 4 (7,m)<0 if 7 G \{V2 - 1)/V2, 1], m = 3,4,5,6. 
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In order to prove (jfi.lHj) and ()6.14j) for each n > 6 we study the sign 
of partial derivative dip^/dn. Set 

(6.19) 6(7, n) = n 2 ~ 2 \2n + l) 2 • ^( 7 ,n); 

on 

then 

(6.20) 6(7 ,„) = -|( 1 + -L)" 2 + 7 ( 1 + _L s " 

7 ^ ' ^ -^n-ir+iri+irVn-^" 1 





V 








r 




n . 




u 




i 






n 



4 \ n. 7 4 \ n y 2 \ n y \ 2n 

11 ^^ivA^v'-^ri- 1 " 



n J \ 2n y 4 \ n 

The power series expansion of 6(7, n) about n = 00 is 



(6.21) 6(7,n) = ^6 & ( 7 )(l/r 

fc=2 

where 

(a oo\ ^ / \ 5 - 22 7 + 18 7 2 -4 7 3 
( 6 - 22 ) 62(7) = 3 1 



(6.23) 



63(7) 



-10 + 25 7 - 14 7 2 + 2 7 ; 



By (I6.20|) . estimating from above |6fe( 7 )|, we obtain 

7 7 1 — 7 



(6.24) 



1 
2 



2 fe 



27I+; 1 



4 2 V2 fc 



+ 7 + 



(k + l)+ 



7 7 



where each term comes from the expansion of the corresponding term 
in IjOty . 

For example, consider 
(6.25) 

7 / 1 \ -1 



1 - 



71 



1 

2ra 



i=i 



j=0 



71 



Since < 7 < 1 we have 
(6.26) '^ 2| 



17 



-(»-!)! < 7_ 
i — 1 ~~ i 



30 PLAMEN DJAKOV AND BORIS MITYAGIN 

Thus the absolute value of the coefficient of {l/n) k in (|fi.25j) does not 
exceed 

1 7 1 7 1 7 1 7 
1_ 1 1_ 1 1_ 1 1 1_ 1 < 

1 7 / 1 1 1 \ 1 

2^ + 2 + 2^ + 2^ + ' ' ' + 1 J ~ ^ + 7 ' 

The inequality (|fi.24j) may be written as 

(6.27) i, i ( 7 )<?.^ + i + i±2 + 2 7 + i^(* + l) + ^. 
Since 

! /I ! A \ / 



(1 -x) 5 



fc=4 fc=4 

we obtain, by f)6.27|) . that 

oo 1 

(6.28) Y,\h(l)\n- k < M^n) ■ ^, 

fc=4 

where 
(6.29) 

w/ , 3(10n-4) 1 +7 337 1~7 5n - 4 

M (7, ™) = 777777 777 + 77777 ^77 + 7777 ,x + 



16(2n-l) 2 8(2n-l) 16(n - 1) 4 (n-1) 2 ' 
Thus we have 

(6.30) n6 2 ( 7 )+6 3 (7)-M(7,n) <n 3 6( 7 ,n) < nb 2 (j)+b 3 (-f)+M^, n). 

On the other hand, 

86 2 ( 7 ) = (5 - 2 7 )( 7 - (1 - l/v / 2))( 7 - (1 + 1/V2)), 

and therefore, 
(6.31) 

62(7) > for 7 G [0, 1 - 1/V2), 62(7) < for 7 G (1 - l/v 7 ^, 1]. 

One can easily see, for each fixed 7 G [0, 1], that Af (7, n) is a decreasing 
function of n. This fact leads, in view of ()6.30|) and (|6.31|) . to the 
following inequalities: 
(6.32) 

< 663(7) + 63(7) -M( 7 , 6) < n 3 6( 7 ,n), 7 G [0,0.19), n > 6 

and 
(6.33) 

n 3 b{^,n) < 66 2 (7)+& 3 (7)+M(7,6) < 0, 7 G (l-l/y/2, 1], n > 6. 
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(We checked the left inequality in (jfi.32j) and the right inequality in 
(J6.33j) numerically by graphing the corresponding functions of 7.) 

In view of (fOTH) and flEEJ, dfi/dn(-f,n) > if 7 G [0,0.19] and 
n > 6, so ^(7, n) increases with n. Therefore, for each 7 G [0, 0.19] and 
n > 6, we obtain by (J6.17|) that 5(7, n) > 5(7,6) > 0, which proves 

dsn. 

In a similar way (|6.33|) implies that (p{p(,n) decreases with n if 
7 G [1 — l/y/2, 1] and n > 6. Thus, in view of (|6.18|) . we obtain that 
5(7, n) < 5(7,6) < 0, for 7 G [1 — l/v2, 1] and n > 6, which proves 
()6.14|) . This completes the proof of Lemma El 

□ 

7. Conclusion; comments and questions 

1. So far in our analysis we focused on the tri-diagonal matrices 
given by (|2.1|) and ()2.2|) . or (jl.3j) with a < 1, or even with a < 1/2. 
The Whittaker-Hill matrices (IH)j) satisfy (jTU) and (j2~2~|) with a = 
1, M = 4 + t. Proposition |U tells us that the eigenvalues E n (z), n G N, 
are analytic functions in the disk A = {\z\ < 1/(8M)}, and nothing 
more. But these matrices come from the differential operator 

(7.1) Ay = -y" + q(x)y, 
considered with 

(7.2) q{x) = a cos 2a; + &cos4r, a = — Azt, b=—2z 2 . 

Let q(x) be a real analytic periodic function of period it. Of course, 
then q extends analytically in a neighborhood of / = [0, 7r], say, in 

(7.3) G £ = {w = x+iy : — s < x < tt+e, — e < y < e}, 3e > 0. 

In other words, q is in the Banach space A(G £ ) of all functions that are 
continuous in G e and analytic in its interior, with the norm 

11/11 =max{|/H|: w G G £ }. 

Consider the boundary conditions 

Per+: y(0)=y(ir), y'(0) = y'(n), 

Per-: y{0) = -y(n), y'(0) = -y'(n), 

Dir : y(0) = y(vr) = 0. 

To be certain, let us talk only about the periodic boundary condi- 
tions Per + , and let us consider the (invariant) subspace of even func- 
tions. Then the operator ()7.1|) has eigenvalue functions E n (z), E n (0) = 
(2n) 2 , n = 1,2,.... 
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H. Volkmer jSHj proved that if q is a real analytic function, then 
E n (z) is well defined as an analytic function in the disk 

A n = {z : \z\ < R n }, R n = an 2 , a > 0. 

Careful analysis of the proof in shows that a stronger quantitative 
statement holds. 

Proposition 18. If 

(7.4) q e A(G £ ), e > 0, 

then the eigenvalues E n (q) of the operator \7. l\j are well defined if q is 
real-valued on [0, tt] and small by norm. Moreover, for each n, E n (q) 
can be extended as an analytic function of q in the ball 

B(R n ) = {qe A(G £ ) : ||g|| < R n }, 

with R n = an 2 , a = a(e) > 0. 

As soon as we have this Proposition, we can consider the potentials 
(17. 2j) as elements of A{G e ), with, say, e — 1/4. Then 

(7.5) ||4^cos2x + 2z 2 cos4s|| < A\zt\e 2 '^ + 2|^| 2 e 4 ^ < 7 (\zt\ + \z\ 2 ) , 
and therefore, if 

(7.6) \tz\ + \z\ 2 <^n 2 , 
then 

(7.7) e n (z) = E n (q), q = Azt cos2x + 2z 2 cos4x, 
is an analytic function of z. Choose 

(7.8) R n = n(l + 4\t\/a)- 1 ; 
then 

(7.9) zeA n = {z:\z\< R n } => (JTSD, 

and therefore, the function e n (z) is analytic in the disk A n . 

We explained the following statement (which is stronger than its 
analogue coming from Proposition 0]). 

Proposition 19. Under the conditions M.b]) the spectrum of the op- 
erator M.l\) is discrete. The function e n (z) G Q7. 7| ) is analytic in 
A n e (T^l, and 

e n (0) = n 2 , \e n (z) — n 2 \ < n if z E A n . 
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2. Of course, the claim of Proposition EE with R n = an/ (a + 4|£|), 
is stronger than Proposition 0] with R n = 1/8. This example, together 
with Remark fl2l supports our belief that, for matrices (L, B) G (|2.1|) . 
Proposition 0] can be significantly improved, so that to give analyticity 
of E n (z) G flUD, (Q in the disk A n G (Q with 

(7.10) = 6n 2 " Q , 36 = 6(a) > 0. 

If a = this is true, but again it comes from H. Volkmer's result |3*7|l3*n] 
for the Mathieu differential operator which is unitary equivalent to the 
matrices (fOjl . (PJ) . 

However, even in this case, no approach to the proof of this statement 
is known in the framework of matrix analysis. 

3. Of course, if the Taylor expansion 

oo 

E n (z) =n 2 + J2 a 2k(n)z 2k 

k=l 

is known, then one may find the radius of convergence of E n (z) as 

r n = I limsup \a 2 k{n)\ l 
V k 

Proposition [7| gives that 

\a2k(n)\ < 8kn ■ 

for (L, B) G (EU) + (Q. However, if B G i.e., 
h = c k = k a } < a < 2, 

we believe that 

(7.11) |a 2fe (n)|<n^f , n = 1, 2, . . . , T > 0, A > 0, 

Of course, fl7~TlT) would imply fl7~TUl) . 

4. Maybe, the representation (|3~Hjl . (I3~2^|) . (jB~2^|) of Propositions □ 
and |H1 could be used in an attempt to get (|7.11|) . But let us make a 
couple of elementary remarks to Propositions [7| and |S1 

Remark 20. It was observed in ()3.8|) . on the basis of the representation 
(|3.6I) and ()3.23|) . f)3.24|) . that a k {n) = for odd k. This follows also from 
the equality 

(7.12) Sp(L + zB) = Sp(L-zB), z G C, 
because \7.1£\ ) implies that all E n {z) are even functions. 





34 



PLAMEN DJAKOV AND BORIS MITYAGIN 



(In particular, this implies that in formulas like 1)4.58)1 and (|4.59|) the 
coefficients Pk(z) should be even functions. In [5], however, formula 
(8) in Theorem 2.1 has P\(z) = (z 3 — 4z)/16, so one can conclude that 
this is not correct even without knowing the correct formula.) 

To get ()7.12|) . consider the unitary operator U defined by 

(7.13) Uej = {-l) j ej, 1 < j < oo, U 2 = 1. 

Then for each matrix A = [A(i,j)} the operator A = U~ l AU = UAU 
has a matrix A(i,j) = (— j). In particular, U~ l (L + zB)U = 
L — zB, i.e., the operators L + zB and L — zB are similar, and therefore, 
()7.12j) holds. Of course, this implies that E n (z) are even functions. 

Remark 21. By Propositions[7|and|Sl the integrals that appear in 
and (|3.24|) vanish if \j — n\ > k. But they vanish even if \j — n\ > k/2. 

After Remark |2*U1 we can talk only about even k, say k = 2m. Let us 
focus on (|3.24j) . i.e., on the integrals 

(7.14) I(n;j,k)= [ \(R° x (BR° x ) k e 3 , e 3 )d\ 

Jh n 

where h n = {A G C : A = n 2 + n + it, t G E}. 

The integrand in ()7.14j) is a linear combination of rational functions 
like ()3.26|) with coefficients depending on B, where each rational func- 
tion corresponds to a walk (jo,ji, ■ ■ ■ ,jk) from j to j on the integer 
grid Z, with steps ±1. Indeed, when the operator R°(BR°) k acts on ej, 
then (since R°e u = (1/(A — v 2 ))e v while Be u is a linear combination of 
e I/ _ 1 and e„ + i) we get a linear combination of 2 k vectors, each of them 
coming from some walk (j , ji, ■ ■ ■ ,jk) as e JO — > — > ■ ■ ■ — > ej k . Since 
(ej k , ej) 7^ only for jk = j, we consider further only walks from j to 
j- 

Moreover, the argument used to prove the point (iii) in the proof 
of Proposition |H1 shows that the rational function Q of (j3.26|) yields a 
non-zero integral over the line h n only if it has poles both on the left 
and on the right of h n , and its poles j 2 come from the vertexes of the 
corresponding walk (jo,ji, . . . , jk) - In other words, if j < n (respectively 
j > n) then the corresponding walk j = ■ ■ ■ ,jk =3 should pass 
through n + 1 (respectively n). 

Take now any j such that \j — n\ > k/2. If j < n (respectively 
j > n), then there is no fc-step walk from j to j passing through n + 1 
(respectively n) because the steps are equal to ±1. Thus each of the 
integrals ()7.14|) vanishes if \j — k\ > k/2. 
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5. We consider a > in (jl.3|) and elsewhere to have unbounded or 
non-compact operators B. Of course, Theorems Q and 121 remain valid 
for a < as well. But then a simpler proof can be given because for 
a < the restriction a < 1 — 2/ k holds with k = 2. In particular, by 
(16. 5|) . i.e., by Proposition [THl we have 

S^(M)= £i4 2) (o) = o 

for any equivalence class of the Spectral Riemann Surface of the pair 
(L,B) e (J2U) + (EH, a < 0. 

Of course, it is easier to study the sign of 0-2(0:, n) than the sign of 
0.4(0:, n) (compare to LemmaH^). By ()3.32|) we have that 



/1X b\C\ . s b n ^\c n _\ b n c n 

a 2 (l) = — , a 2 (n 



3 v ' 2n-l 2n + l 

If (b n ) and (c n ) are decreasing sequences of positive numbers, then 

a 2 (l) < 0, a 2 (n) > 0, n > 2, 

and we can use the same argument as before (see the proof of Theorem 
0I) to conclude that the corresponding Spectral Riemann Surface is 
irreducible. So, we proved the following analogue of Theorem El 

Proposition 22. Suppose that \2. 1)) and JjQ| ) hold with monotone 
decreasing sequences b = (b n ) and c = (c n ), and with a < 0. Then the 
corresponding Spectral Riemann Surface is irreducible. 

We have to admit that with all variety of pairs (L, B) for which 
we have proved the SRS's irreducibility, we know no nontrivial (i.e., 
beside the case where some entries bk or Ck vanish, or diagonal entries 
are multiple) example of a pair (L, B) with a reducible SRS. 

6. From a < we can go to another direction, i.e., consider a G 
(1/2,1). The estimate (|7.11|) is our conjecture, but even now we can 
claim the following amendment to Theorem ^ 

Proposition 23. Under the assumptions \2.1\) and \2.2}) . if < a < 

9/10, then the regularized trace 

(7.15) tn(z) = f; (E n {z) -n 2 - ^E'MzA 

n=l ^ ' 

is well defined as an entire function of z, and 

(7.16) tn{z) = 0. 
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The proof is based on 1)4.41)1 - 1)4.44)1 and the estimates given by 
Lemma It goes along the same lines as Definition of regularized 
trace in Section 5.4 and the proof of Theorem [TJ see ()5.20)) to ()5.22j) . 
We omit the details. 

Of course, one can introduce the higher order regularized traces 



and study for which a this expression is well defined as an entire func- 
tion. 

It is important to mention that many interesting examples of eval- 
uation of a regularized trace can be found in the recent papers |HJ 
U EE I2D 122 123 121 although there the operators L and B are usually 
self-adjoint and z is real. Let us notice that, in our Theorem Q the 
first line of (jl.7|) . a < 1/2, can be interpreted as an example to Thm 1 
in [22] • Then, the second line of (jl.7j) shows that the restrictions on 5 
and lj in [22 a could not be weakened. 
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